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We present a model for host-parasite dynamics which incorporates both vertical and horizontal 
transmission as well as spatial structure. Our model consists of stacked contact processes (CP), 
where the dynamics of the host is a simple CP on a lattice while the dynamics of the parasite is a 
secondary CP which sits on top of the host-occupied sites. In the simplest case, where infection does 
not incur any cost, we uncover a novel effect: a nonmonotonic dependence of parasite prevalence on 
host turnover. Inspired by natural examples of hyperparasitism, we extend our model to multiple 
levels of parasites and identify a transition between the maintenance of a finite and infinite number 
of levels, which we conjecture is connected to a roughening transition in models of surface-growth. 
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I. Introduction 

The need to understand and control the dynamics of 
infections has motivated the development of a variety of 
statistical mechanical models. One of the most important 
of these is the contact process (CP) PHI], which describes 
the dynamics of an infection in a spatially structured pop- 
ulation. In the contact process, each site on a lattice rep- 
resents a host organism which can be infected or suscepti- 
ble (uninfected) . Infection is transmitted to neighbouring 
susceptible host sites at rate b, and infected hosts recover 
(i.e. become susceptible) at rate d. The contact process 
provides a baseline model for many problems in ecology 
and epidemiology. It is also of fundamental importance 
in statistical physics. This is because it exhibits a non- 
equilibrium phase transition between an infected and a 
non-infected phase at a critical value of A = 6/d 

The standard CP model for an infected population as- 
sumes that the host population is of fixed size, without 
turnover, i.e. host births and deaths. This is valid if 
the timescale on which the infection is gained and lost is 
much shorter than the lifespan of an individual. Some 
infections, however, are carried by individuals for long 
times, and may be transmitted to offspring upon repro- 
duction (vertical transmission) , as well as being transmit- 
ted horizontally — i.e. upon physical contact between in- 
dividuals. Examples include plasmids carried by bacteria 
, some microsporidian parasites of insects Q , maternal 
transmission of parasitic worms (helminths) in humans 
and other animals Q , and animal viruses including HIV 
and several hepatitis viruses. Here, we investigate the 
interplay between vertical and horizontal transmission in 
such populations. Other authors have determined the 
conditions for parasite persistence in mean-field models 
that lack spatial structure, in which the parasite may 
affect host fitness [8l4l0j. and have suggested that these 
conditions may be affected by spatial structure [ll[ . Here, 
we show that spatial structure can produce a qualitatively 
new effect: a coupling between the dynamics of the in- 
fection and of the underlying host population, even when 
the infection does not affect the fitness of the host. 

We present a two level stacked contact process, in which 



the host population is represented by a CP on a lattice, 
and the parasite population is represented by a second 
CP which sits on top of the host CP. This model incor- 
porates both vertical and horizontal transmission — if 
an infected host reproduces or dies, the parasite is re- 
produced or dies with it (vertical transmission), and a 
parasite-infected site can infect a neighbouring site if it 
is occupied by an uninfected host (horizontal transmis- 
sion). We characterize the conditions for parasite per- 
sistence in the form of a phase diagram and discover an 
interesting phenomenon: although the steady-state prop- 
erties of the host population depend only on the ratio of 
its birth and death rates, A = b/d, the prevalence of the 
parasite can depend non-monotonically on the host pop- 
ulation's turnover rate. This phenomenon has its origins 
in the fluctuations of spatial clusters of host individuals, 
and cannot straightforwardly be captured by a mean-field 
theory. 

Parasitic infections are not always limited to two levels. 
Hyperparasitism, in which an organism carrying a pri- 
mary parasite is susceptible to a secondary parasite [12j . 
can be harnessed as a biocontrol mechanism — exam- 
ples include viral infections of the fungus Cryphonectria 
parasitica that causes chestnut blight [13J, and cytoplas- 
mic RNA elements that infect the fungus causing Dutch 
elm disease Indeed, this parasites on parasites sce- 



nario was envisaged more than 250 years ago by Jonathan 
Swift, who wrote: "So, naturalists observe, a flea, has 
smaller fleas that on him prey, and these have smaller still 
to bite 'em, and so proceed ad infinitum" [l5|. Inspired by 
this scenario, we extend our model to a multilevel stacked 
contact process, in which individuals carrying a primary 
infection may be susceptible to secondary infections, and 
those carrying the secondary infection may be susceptible 
to tertiary infections, etc. This raises a number of ques- 
tions: How does the dynamics of one parasite level couple 
to the next? How many levels of parasites are sustain- 
able in a population — does the hierarchy really proceed 
11 ad infinitum"! We show that in our stacked contact 
process, there is a well-defined transition between main- 
tenance of an infinite hierarchy of levels of parasites, and 
limitation to a finite number of levels. The transition 
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between these two regimes appears to be connected to a 
rou ghe ning transition in certain surface growth processes 
[lg-[l8| . This work presents new challenges to our under- 
standing of contact processes, with potential implications 
for the dynamics of long-timescale infections. 

II. Two-level stacked contact process 

We begin by considering the dynamics of a host and 
a single parasite which can be transmitted vertically or 
horizontally — a two-level stacked CP — on a 2D square 
lattice. Lattice sites can be either empty or occupied by 
a host organism. Occupied lattice sites become empty at 
rate do, due to death of hosts, and host organisms at- 
tempt to replicate into neighbouring empty sites at rate 
60 (the attempt rate in any direction is bo/z, z being the 
coordination number of the lattice) . Host organisms may 
also carry a parasite, which we assume to be neutral in 
the sense that organisms with and without the parasite 
have identical birth and death rates. Infected individu- 
als pass the parasite to uninfected neighbours at rate 61, 
and the parasite is lost at a constant rate d± . We assume 
that the offspring of an infected organism is also infected 
(vertical transmission), and that when an infected host 
dies, the parasite dies with it. Thus the infection dynam- 
ics consists of a secondary CP, with parameters b\ and 
di, which sits on top of the host's birth-death CP. Since 
the host population does not occupy the whole lattice, 
the secondary CP does not take place on a regular lat- 
tice but rather on the irregular network of occupied sites 
which changes stochastically due to the dynamics of the 
host population. The processes that occur in our model 
are illustrated in Figure [TJ We simulate this model us- 
ing a stochastic kinetic Monte Carlo algorithm, in which 
all events occur as Poisson processes [lj|, on a 100 x 100 
lattice with periodic boundary conditions. Averages are 
calculated after an initial transient to allow the steady 
state to be reached. 

III. Host dynamics influences parasite persistence 

In this model, the host population undergoes a simple 
CP, whose stationary properties are fully determined by 
the single parameter Ao = bo /do [IH3|- If Ao is less than 
a critical value, A cr ;t, the only steady state is an empty 
lattice and the population rapidly becomes extinct. If 
Ao > A cr it, a non-zero population can be maintained 
for long times, before eventually becoming extinct due 
to a rare fluctuation. The transition between these two 
regimes is second order; the values of A cr i t and the associ- 
ated critical exponents have been characterised in detail 

3]. Defining r, — or 1 if site i is empty or occupied 
respectively, the dynamics of the host population obeys 



d7 to) = \ E 5 ° to ^ ~ r *)) ~~ d ° to) , 



(i) 



j 



where the sum is over the z neighbours of lattice site 
i and the angle brackets denote averages over multiple 
realizations of the dynamics. The steady-state solution 
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Figure 1. The events that constitute the 2- level stacked CP. 
Red and black circles represent host organisms and parasites 
respectively. A site can be empty, occupied by a susceptible 
host or occupied by an infected host. The left panels represent 
reproduction of a susceptible or infected host (top and middle, 
rate bo), and transmission of the parasite (bottom, rate 61). 
The right panels represent death of a susceptible or infected 
host (top and middle, rate do), and loss of the parasite (bot- 
tom, rate di). z is the number of nearest neighbour sites. The 
dotted squares indicate the level of the CP at which the event 
happens. 



of Eq. (QJ satisfies (r,) = (Ao/(Ao — 1)) (^iTj)- Ignoring 
correlations between neighbouring sites by assuming that 
(tsTj) ~ (rj) (tj), we arrive at the mean- field result (r) = 
1 - (1/Ao). 

Considering now the dynamics of the parasite, we define 
cr, = 1 if site i contains an infected host and <7j = 
otherwise (i.e. if site i is either unoccupied or contains a 
susceptible host). Oi obeys 
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The first term on the r.h.s. of Eq. @ corresponds to ver- 
tical transmission: an empty site is filled by replication of 
an infected host. The second term represents horizontal 
transmission: a susceptible host, denoted by (r, — <Zj), is 
infected by a parasite-carrying neighbour. The final two 
terms correspond to death of an infected host and loss of 
the parasite. Eq. @ can be rewritten as 

^ to) = \ E t 5 ° to) + fa - M (TiCTj) - 61 fa<7j)] 
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Comparing this with Eq. (|T|), we see that the cross- 
correlation (tjCTj) perturbs the parasite dynamics from 
that of a standard CP. Interestingly, however, if b\ = bo 
(i.e. if the rates of horizontal and vertical transmission 
are equal), this term vanishes and the form of Eq. ^ 
becomes that of a standard CP with parameter A e g = 
AoAi/(Ao + Ai) = b/^2 k dk- The parasite dynamics also 
has CP-like behaviour in the mean- field limit, where we 
neglect spatial correlations: setting {naj) « (t») (aj) and 
fatTj) w (a^ ((Tj) in Eq. ((3|), we obtain the steady-state 
solution (cr) = 1 - (1/Ao) - (1/Ai) = 1 - (1/A cff ). 
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Figure 2. Phase diagram of the two- level stacked CP. Sym- 
bols: simulation data showing boundary of parasite persis- 
tence, for di = 1 and do = 1 (circles) or do = 20 (squares). 
Dashed line: mean-field prediction Ai = Ao/(Ao — 1). 



Figure shows the phase diagram for our model, as a 
function of Ao and Ai. Three steady state scenarios are 
possible: (i) the host population is extinct, (ii) the host 
population is finite but the parasite is extinct, and (iii) 
the parasite persists within a finite host population. Since 
the dynamics of the host population is a standard CP, the 
host population is extinct (scenario (i)) if Ao < A cr it- For 
Ao > Adit, the host population persists, with or without 
the parasite. The condition for parasite persistence (i.e. 
the boundary separating scenarios (ii) and (iii)) is pre- 
dicted by the mean field theory to be Ai = Ao/(Ao — 1) 
(obtained by setting (er) = 0) - shown by the dashed 
line in Figure [2j As one might expect, as Ao decreases, 
the density of the host population decreases, and a higher 
rate of horizontal transmission (Ai) is needed to maintain 
the parasite. 

The symbols in Figure [5] show the boundaries between 
parasite persistence and loss obtained from our kinetic 
Monte Carlo simulations, for two values of do, the host 
death rate. The fact that these are shifted upwards and to 
the right of the mean-field prediction shows that spatial 
correlations make it harder to maintain the parasite. It 
is also interesting that we obtain different results for the 
two different values of do. This shows that the parameters 
Ao and Ai do not fully determine the phase behaviour of 
the system. For fixed Ao = bo /do, a higher turnover rate 
of the host population (i.e. higher do and &o) apparently 
makes it harder to maintain the parasite. 

In fact the situation is even more complex. Figure I3al 
shows a more comprehensive investigation of the steady- 
state parasite density (a) as a function of the host 
turnover rate. In these simulations, we vary bo and do 
keeping Ao = b /do fixed. The steady-state density of the 
host population is constant, since it depends only on Ao- 
Remarkably, for some parameter combinations, the den- 
sity of the parasite actually depends non-monotonically 
on the host turnover rate. If the host dynamics is slow, in- 
creasing the turnover rate increases the parasite density, 



while if the host dynamics is fast, the parasite density 
decreases with the turnover rate. 

This dependence of the parasite density on the host 
turnover rate is absent in the mean-field theory and must 
therefore be a consequence of spatial correlations. To 
investigate this, we measure in our simulations the spa- 
tial correlations of the parasite density via the function 
C'i.j = [( <x i cr j)/( <T )] ~~ ( a )- The first term is the condi- 
tional probability of finding a parasite at site j given that 
there is one at site i, while the second term ensures that 
Cij — > as \i — j\ — > oo (where \i — j\ denotes the dis- 
tance between sites i and j). We then fit our data to the 
functional form Cjj ~ e - ^ -3 '''^ to extract the spatial cor- 
relation length £, which can be thought of as the parasite 
cluster size. This is shown in Figure I3bl Strikingly, the 
correlation length is minimal at the turnover rate where 
the parasite density is maximal. This makes intuitive 
sense: since horizontal transmission requires contact be- 
tween infected and uninfected hosts, minimal spatial clus- 
tering of the parasite population results in maximal par- 
asite density. 

Why should the turnover rate of the host affect the spa- 
tial clustering of the parasite? At fast host turnover rates, 
we observe in our simulations that local patches of host 
organisms rapidly grow from a single seed and vanish 
by stochastic extinction, on a faster timescale than that 
of horizontal transmission. Because offspring are always 
of the same type as their parents, organisms in a sin- 
gle patch are either all infected or all uninfected. Thus 
the infected population is highly clustered, leading to low 
rates of horizontal transmission and low parasite density. 
By contrast, when the host population turnover is slow, 
the parasite dynamics constitutes a CP on the effectively 
frozen, disordered network of lattice sites that are occu- 
pied by host organisms. This network contains clusters of 
sites that may be disconnected, or poorly connected, from 
the rest of the host population; this creates the possibil- 
ity of local, stochastic parasite extinctions. This is most 
obviously the case in the low Aq regime where the host 
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Figure 3. (a) Parasite density (a) as a function of host 
turnover rate do for fixed Ao (and hence fixed host den- 
sity). Curves top to bottom correspond to (Ao,6i,di) = 
(5, 4, 1),(2.5, 15, 1),(2, 8, 1) and (1.8, 15, 1). (b) Parasite den- 
sity (a) (triangles, left axis) and spatial correlation length £ 
(squares, right axis) for (Ao, 6i, di) = (1.8, 15, 1). 
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Figure 4. Possible transitions for the multi-level CP, consid- 
ering a site with label m = 3. The symbols are as in Fig. [TJ 
the green circles represent the secondary parasite. Birth pro- 
cesses are shown on the left. The organism in question can 
reproduce into a neighbouring empty site (top left), transmit 
its primary parasite to a susceptible neighbour (middle left 
- note the secondary parasite is also transmitted), or trans- 
mit its secondary parasite to a susceptible neighbour (bottom 
left). Death processes are shown on the right. These consist 
of death of the host (top right), loss of the primary parasite 
(middle right) or loss of the secondary parasite (bottom right). 



density is low, but it is still true for higher host densi- 
ties. In this regime, an increase in the host turnover rate 
serves as a mixing mechanism, by more homogeneously 
distributing the host population and thereby providing 
greater opportunity for the parasite to spread. 
Thus the nonmonotonic dependence of the parasite den- 
sity on the host turnover rate arises from a coupling of 
the spatial fluctuations of the underlying host CP to the 
parasite dynamics, and can be viewed as a competition 
between the mixing effect of the host birth-death pro- 
cess at low turnover rates and the population segregation 
arising from parent-offspring clustering at high turnover 
rates. 

IV. Multi-level stacked contact processes 

Motivated by natural examples of hyperparasitism, we 
now investigate a multi-level stacked CP, where individ- 
uals carrying a primary infection are susceptible to a sec- 
ondary infection, and those with the secondary infection 
are susceptible to tertiary infections, and so on. In par- 
ticular, we are interested in whether or not such a system 
can sustain an infinite hierarchy of parasite levels ("ad 
infinitum" ) . 

In the M-level stacked CP, a site is labelled if it is 
empty, 1 if it contains an uninfected host, and m = 
2, 3, . . . , M if it contains a host infected with m — 1 lev- 
els of parasites. For example, a site labelled 2 contains 
a host infected with only a primary infection (parasite) 
while a site labelled 3 contains a host which is infected 
with primary and secondary infections (a parasite and a 
hyperparasite) . Note that the presence of an infection at 
level m implies the presence of lower- level infections (e.g. 
one cannot have a secondary infection without a primary 
infection) . Figure @] illustrates the dynamical processes 
that can occur in this model. An empty site may be oc- 



cupied at rate bo/z by reproduction of a neighbouring site 
of any label (top left), or a site labelled m may be infected 
at rate b m /z by higher-level parasites from a neighbour- 
ing site with label n > m; its label is then promoted to 
n (middle and bottom left). The host organism of a level 
m site can die, at rate do (top right), or experience loss 
of one of its parasites, at rate d n (where n < m relates 
to the level of parasite that is lost). When the latter 
happens, all higher parasites are also lost and the site is 
demoted to level n (middle and bottom right panels in 
Fig. 0]) . The parameters of the model are the number of 
levels M, and the set of level-dependent birth and death 
rates {6 , ■ . .,b M -i;d , ■ . . ,d M -i}- 

We hrst consider the special case where the birth rates 
at all levels are equal; b n = b V n. We divide lattice sites 
into two sets: sites that have labels greater than or equal 
to m (here denoted M+), and sites that have labels less 
than m (here denoted An M- site can become 

M.+ by infection by a neighbouring ftA + site: this occurs 
at rate b/z regardless of the level of the M.+ site. An 
M.+ site can become M.- by death of the host or by loss 
of a parasite at levels n = 2 . . . m. This occurs at total 
rate d = X^flo 1 ^ where the sum is over levels up to 
m—1. Thus the dynamics of the M. + density is exactly 
equivalent to that of a standard CP with parameter A e g = 

6/ zXo 1 d e- 

This observation provides us with important insight into 
the number of levels that are sustainable in a stacked 
CP. For the standard (one-level) CP, a finite density of 
occupied lattice sites can only be sustained for A > A cr ;t, 
where A cr it ~ 1.649 for a 2D square lattice ||. In the 
stacked CP with equal birth rates we therefore expect the 
density of M.+ sites, those with label n > m, to be non- 
zero only if A c ff = bo/ Y^eLo 1 < A cr it ■ Figure [5] shows 
simulation results for the average density of lattice sites 
with label n > m, as a function of m, for a stacked CP 
with equal birth rates b n = b and death rates d n = d = 1, 
for several values of b. By mapping onto a standard CP 
with A c ff = b/(md), we predict that all levels n > m* = 
\b/ (A cr itG?)l have zero density and are thus unsustainable. 
Our simulation results, shown in Figure [5l bear this out: 
the system indeed only sustains a finite number m* of 
parasite levels. 

Are there any circumstances where a stacked CP can 
be sustained 11 ad infinitum" — i.e. for an infinite number 
of levels? This is indeed possible if either the birth rate 
b n increases, or the death rate d n decreases, sufficiently 
strongly with n. We first suppose that the birth rate 
b n = b is constant but the death rate decreases by a 
factor / at successive levels: d n = fd n -\ = f n do where 
< / < 1. In this case, the average density of sites 
with label n > m is given by that of a standard CP with 
Acff = bj ^Xo 1 <k = 6(1-/)/ - f m )}- To sustain an 
infinite number of levels, we require that AJ^ 00 = 6(1 — 
f)/do > A cr it, which implies that / < 1— (doA C rit/6). Our 
simulation results, Fig. I6al show that, with do = 1 and 
6 = 5, the system indeed sustains an infinite number of 
levels for values of / above the critical value of / ~ 0.67. 
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Figure 5. Density of sites with label m or higher, as a 
function of level m, for a stacked CP with equal birth rates 
b n — b and death rates d n = 1. Left to right correspond 
to simulations with b — 2,4,6,10 and 15 respectively. Arrows 
indicate m* = |~&/(Acritd)] , the boundary for sustainability as 
predicted by mapping to the standard CP. 




Figure 6. (a) Density of sites with n > m, as a function of 
level m, in the stacked CP for decreasing death rates d n — 
f"do, with do = 1, b n = 5 (solid lines) and for increasing birth 
rates b n = f~ n bo, with bo = 5, d n = 1 (dashed lines). In 
the former case, the exact CP predicts a nonzero density as 
m — ¥ oo only for / ^ 0.67. (b) Density of the m th level in the 
limit m — > oo with b n = 5, d n = /"do, / = 0.75 as a function of 
host death rate do from simulations (circles) and as predicted 
by a standard CP (triangles) with A"^°° = 6(1 - f)/d . The 
simulation data is obtained by measuring the plateau values 
in plots like panel (a), for m values up to 20. The mapping 
to a standard CP predicts that only a finite number of levels 
are sustainable for do > dg « 0.758. 



Interestingly, this analysis also shows that, for any given 
values of b and /, there exists a critical value of the host 
death rate do which separates regimes where the stacked 
CP can and cannot sustain an infinite number of levels. 
This is given by = 6(1 — /)/A cr i t . Figure IFbl shows, 
as a function of do, the predicted density of sites with 
label n > m, as m — »• oo: this is given by the density of 
a standard CP with A c ff = b(l — f)/do. The density as 
m — > oo indeed falls to zero at do = d^. 

Finally, we explore briefly the case where, instead of 
varying the death rate between levels, we instead increase 
the birth rate by a factor 1/ f {i.e. set b„ = b n -\/f). 
In this case the dynamics no longer maps onto that of 
a standard CP - but nevertheless, our simulations show 
qualitatively similar results (dashed lines in Figure I6ap . 



The transition to an infinite number of levels of para- 
sites that we observe in this model appears to be related 
to a transition that occurs in models for the growth of in- 
terfaces by deposition of particles on surfaces [Tojfisj . In 
these models, the interface is modelled as a lattice, with 
a given height at each lattice site. According to the bal- 
listic deposition rule for surface growth [2(| , particles fall 
onto the lattice from above and only stick if the neigh- 
bouring lattice site already contains a particle. In some 
models, particles can also desorb from the surface fl8j |. 
In our parasite model, we can think of the label m of a 
given lattice site as corresponding to the local height of 
the interface. The transmission of higher-level parasites 
to a site with label m, from a neighbour with n > m 
then corresponds to ballistic deposition, while death of 
a host and loss of parasites loosely correspond to des- 
orption of particles. This apparent mapping is intriguing 
because these models for surface growth show a rough- 
ening transition [i~fi - tl8| : when the deposition rate is low, 
the interface remains smooth, that is, the width of the 
surface layer remains finite in the thermodynamic limit, 
whereas when the deposition rate is high, the surface layer 
grows and roughens over time, that is, arbitrarily large 
differences in surface height can arise in the thermody- 
namic limit. These two cases correspond to finite and 
infinite hierarchies of parasites in our stacked CP model. 
Although the mapping between our model and these sur- 
face growth models is not exact, one expects to see the 
same phenomenology, e.g. in terms of critical exponents, 
for the transitions. 

V. Discussion 

In this paper, we have presented a simple model for 
the dynamics of long-lived infections in spatially struc- 
tured populations. The model is a stacked contact pro- 
cess (CP), in which the host population undergoes a stan- 
dard CP, and a secondary CP representing the parasite 
takes place on the dynamically changing network of lat- 
tice sites occupied by the host. We find that the cou- 
pling between the dynamics of the host population and 
of the infection leads to non-trivial effects, including a 
non-monotonic dependence of the parasite prevalence on 
the host population's turnover rate, for fixed host pop- 
ulation size. This exposes a connection between spatial 
and temporal scales, since host population turnover af- 
fects spatial clustering of infected organisms which in turn 
affects the prevalence of the infection within the popula- 
tion. We have found that an improved mean-field theory 
that takes local correlations into account (the ordinary 
pair approximation (2l| ) fails to reproduce this nonmono- 
tonic behaviour. This suggests that a deep understanding 
of the fluctuations in the host CP is required to predict 
the steady-state properties of the parasite population. It 
will be interesting in future to investigate how the nature 
of the CP phase transition for the infection dynamics is 
affected by the dynamics of the underlying host popula- 
tion. Since disorder can alter critical exponents and lead 
to new dynamics such as Griffiths phases and activated 
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scaling [22], |23| , we envisage that the secondary CP may 
show fundamentally different behaviour to the standard 
CP. 

Inspired by natural examples of hyperparasitism, or 
"parasites on parasites" , we have also investigated the 
properties of a multi-level contact process. For this 
model, we show that the average density of the m-th 
(and above) level of the stacked CP maps onto that of 
a standard, one-level CP, if the birth rates at all levels 
are equal. We find that a phase transition separates two 
very different behaviours of the model: maintenance of 
an infinite hierarchy of levels of parasites, versus collapse 
to zero density at a finite number of levels. In light of 
the link to models of surface growth it will be interest- 
ing to determine more fully the properties of the parasite 
maintenance transition that we observe in the multi-layer 
stacked CP. 

Finally, we note that our model also presents the possi- 
bility for a novel study of disorder in contact processes. 
Classic models for disordered systems usually create dis- 



order by removing a fraction of sites from the dynamics 
24| or by allowing different creation or annihilation rates 
for different sites (25|. This disorder is often quenched 
- i.e. time- invariant. In our model, disorder in the sec- 
ondary CP (and higher CPs if present) arises naturally 
from the underlying dynamics of the primary CP, and can 
be controlled by varying the rates. This naturally arising 
disorder is not fixed but varies in time and space (in a fun- 
damentally different process than the mobile disorder of 
Ref . [Ill ) . Future studies of the effects of this natural dis- 
order on the properties of the secondary CP may uncover 
new principles of CP dynamics on disordered lattices. 
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